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Abstract. Time resolved spectroscopy of the Intermedi- 
ate Polar FO Aqr reveals the presence of multiple peri- 
odicities in the UV range. A strong orbital modulation 
dominates both continuum and emission line flux vari- 
abilities, while line velocity motions are only detected at 
the rotational frequency. A prominent orbital periodicity 
is also observed in coordinated optical photometry, where 
FO Aqr was previously found to be spin dominated. The 
spectral dependence of the main periodicities shows the 
presence of multi-temperature components in FO Aqr and 
for the first time a hot and a cool component in the rota- 
tional modulation. From a comparison with previous UV 
and optical data obtained in 1990, no spectral variations 
in the orbital and rotational variabilities are detected, in- 
dicating no significant changes in the effects of X-ray illu- 
mination but rather a shrinking of the accretion curtain 
accompained by an increase in size of the thickened part of 
the accretion disc. These observations, consistent with the 
recently discovered long term trend in the X-ray pulsation 
amplitudes, independently confirm a change in the accre- 
tion mode in FO Aqr, which switched from a disc- fed into 
a disc-overflow state, likely triggered by mass accretion 
variations. 

Key words: accretion - binaries: close - stars, individual: 
FO Aqr - Ultraviolet: stars X-rays: stars 



1. Introduction 

Intermediate Polars (IPs) are a subclass of magnetic Cat- 
aclysmic Variables which consist of an asynchronously ro- 
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fating (Pspin < Porb) magnetized white dwarf accreting 
from a late type, main sequence, Roche-lobe filling sec- 
ondary star (Patterson 1994; Warner 1995). 

Except for a few systems for which polarized opti- 
cal/IR emission is detected, the white dwarf is believed 
to possess a weak (< 2 MG) magnetic field which domi- 
nates the accretion flow only at a few radii from its surface. 
Within the magnetospheric radius, material is channeled 
towards the magnetic polar regions in an arc-shaped ac- 
cretion curtain (Rosen et al. 1988). At larger distances, 
different accretion patterns can be present: a truncated 
accretion disc (disc- fed systems) , direct accretion from the 
stream onto the magnetosphere (disc-less systems) as well 
as a combination of the two, where the stream material 
overpasses the disc (disc-overflow) (Hellier 1995 and ref- 
erences therein). 

Due to the asynchronous rotation, IPs show a wide 
range of periodicities at the white dwarf spin (cj), the 
orbital (il) and sideband frequencies (Warner 1986; Pat- 
terson 1994, Warner 1995), whose amplitudes can be 
different in different spectral ranges (de Martino 1993). 
FOAqr (112215-086) was known to show strong periodic 
X-ray, optical and IR pulsations at the spin frequency, 
u = l/(Pspin — 20.9 min), and lower amplitude varia- 
tions at the orbital = l/(Poib — 4.85 hr) and beat 
u! — ft — l/(Pboat — 22.5 min) frequencies (de Martino et 
al. 1994, hereafter Paper 1, and references therein; Marsh 
& Duck 1996; Patterson et al. 1998). The dominance of the 
spin pulsation at optical and high X-ray energies (>5 keV) 
can be accounted for by a disc-fed accretion, whose evi- 
dence was provided by a partial eclipse in the optical con- 
tinuum and emission lines (Hellier et al. 1989; Mukai et 
al. 1994; Hellier 1995). FOAqr was also found to possess 
a disc-overflow accretion mode (Hellier 1993), and the re- 
cent evidence of a long term variability in the amplitudes 
of the X-ray pulsations has been interpreted as changes 
in the accretion mode (Beardmore et al. 1998). This kind 
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of variability, only recently recognized, is also observed in 
other IPs like TXCol (Buckley 1996; Norton et al. 1997) 
and BGCMi (dc Martino ct al. 1995). 

The identification of the actual accretion geometry and 
the determination of energy budgets of the primary X-ray 
and secondary reprocessed UV, optical and IR emissions 
reside on multi- wavelength observations. The main modu- 
lations in FO Aqr have been studied in the X-rays and at 
optical/IR wavelengths. The optical spin pulsations, oc- 
curring mostly in phase with the X-ray ones, were found 
to arise from the outer regions of the accretion curtain 
(Paper 1; Welsh & Martell 1996). The orbital modulation 
from UV to IR was instead found to be multicomponent 
and attributed to the X-ray heated azimuthal structure 
of the accretion disc (henceforth bulge) and to the inner 
illuminated face of the secondary star (Paper 1). However 
the temperature of the UV emitting bulge could not be 
constrained due to the lack of a precise quantification of 
the spectral shape of the UV spin modulation. In this work 
we present high temporal resolution spectroscopy acquired 
with HST/FOS which provides the first detection of dif- 
ferent UV periodicities in FOAqr. For a comprehensive 
study, these data arc complemented with low temporal 
resolution lUE spectra along the orbital period. Coordi- 
nated optical photometry extends the study to a wider 
sp(;ctral range, providing the link to investigate the long 
term behaviour of these variabilities. 

2. Observations and data reduction 

The UV and optical campaign on FO Aqr was carried out 
between September and October 1995 with HST, lUE and 
at the South African Astronomical Observatory. The jour- 
nal of the observations is reported in Table 1. 

2.1. The HST data 

HST Faint Object Spectrograph observations of FOAqr 
were performed on September 10, 1995. The observations 
were carried out in the rapid mode during 7 consecutive 
HST orbits. Due to target acquisition procedures the to- 
tal effective on source exposure time was 4.07 h, yielding 
six continuous exposure slots as detailed in Table 1. The 
orbital period was unevenly sampled since it is commen- 
surable with the HST orbit. The G160L grating was oper- 
ated with the blue digicon covering the range 1154-2508 A 
at a resolution of 6.8 A diode"^ and with the 0.86" upper 
square aperture, supposed to be free from the 1500-1560 
A photocathode blemish which is known to affect the cir- 
cular apertures (HST Data Handbook, 1997). A total of 
797 spectra were collected, each with an effective exposure 
time of ~ 18 s. 

The standard routine processing STScI pipeline ap- 
plied to the data at the time of the observations revealed 
the presence of anomalous features when comparing the 
reduced spectra with the lUE data, and in particular in 



Table 1. Journal of observations. 



Instrument 


Date 






UTstart 

hh:mm:ss 


Duration 
s 


HST/FOS 












Slot 1 


1995 Sept 


. 10 


14:11:04 


2332 


Slot 2 


1995 Sept 


. 10 


15:45:11 


2461 


Slot 3 


1995 Sept 


. 10 


17:21:18 


2461 


Slot 4 


1995 Sept 


. 10 


18:57:49 


2461 


Slot 5 


1995 


Sept 


. 10 


20:34:22 


2461 


Slot () 


i!)!).") 


Soi)t 


. iO 


22:10: ") i 


2 1()1 


lUE 












Image # 












o VV r^ODUoy 


1995 


Oct. 


10 


1 A•AA'^ Q 

14.44. lo 


zozu 


o VVr^ODUyU 


1995 


Oct. 




lO.zl.zD 


zozu 


O VV r^ODUyi 


1995 


Oct. 


10 


lo.Ul.l I 


zozu 


o VV Jr^ODUyz 


1995 


Oct. 


1 K 

10 


1 Q.OO.QQ 

ly. zz. oy 


zozu 


o VV i ODuyo 


1995 


Oct. 


10 


ZU.ao.OD 


zozu 


D VV i oouy4 


1995 


Oct. 


10 


oo.n9.07 
ZZ.UZ.Z / 


zozu 




1995 


Oct. 


10 


Zo.zy. 01 


zozu 


SWP56096 


1995 


Oct. 


16 


00:57:16 


2520 


SWP56097 


1995 


Oct. 


16 


02:24:39 


2520 


SWP56098 


1995 


Oct. 


16 


03:55:35 


2520 


LWP31592 


1995 


Oct. 


15 


13:45:09 


2520 


LWP31593 


1995 


Oct. 


15 


15:34:55 


2520 


LWP31594 


1995 


Oct. 


15 


17:17:11 


2520 


LWP31595 


1995 


Oct. 


15 


18:45:29 


1260 


LWP31596 


1995 


Oct. 


15 


20:06:28 


1260 


LWP31597 


1995 


Oct. 


15 


21:25:29 


1260 


LWP31598 


1995 


Oct. 


15 


22:46:08 


2520 


LWP31599 


1995 


Oct. 


16 


00:13:32 


2520 


LWP31600 


1995 


Oct. 


16 


01:41:01 


2520 


LWP31601 


1995 


Oct. 


16 


03:26:04 


1260 


SAAO 












BVRI 


1995 Oct. 


18 


18:37:44 


13758 




1995 Oct. 


19 


19:41:11 


2745 




1995 Oct. 


20 


20:01:30 


7480 


No Filter 


1995 Oct. 


21 


18:04:35 


15585 




1995 Oct. 


23 


18:40:20 


9130 



the regions 1500 1590 A and 1950 2010 A . The data were 
then re-proccsscd using the STSDAS/CALFOS routine 
within IRAF using the latest reference files for sensitiv- 
ity correction and appropriate aperture flatfield provided 
by the ST Scl Spectrograph Group in summer 1998. The 
calibrated spectra appear then free from the above fea- 
tures although the lUE fluxes are on average ~ 1.1 times 
larger than the EOS ones (see Fig. 1). A check against sys- 
tematic effects using spectra of standard stars observed 
with both lUE and FOS only confirms the known lower 
fiux (on average of ^ 7%) of lUE, with respect to that of 
FOS (Gonzalez- Riestra 1998). The residual flux difference 
may be due to the fact that the lUE data were acquired 
a month later than the FOS spectra. 



de Martino et al.: Time Resolved HST and lUE 



3 



The G160L grating provides zero order light covering 
the full range between 1150 and 5500 A, with an effective 
wavelength at 3400 A, which provides useful simultaneous 
broad-band photometry. The flux calibration of the zero- 
order signal (Eraclcous & Horne 1994) updated for errors 
and post-COSTAR sensitivity and aperture throughput 
(HST Data Handbook 1997) has been applied to the sig- 
nal extracted in the zero-order feature of each of the 797 
exposures. 

2.2. The lUE data 

On October 15, 1995 ten lUE SWP (1150-1980 A) and 
ten LWP (1950-3200 A) low resolution (~ 6 A) spectra 
were acqiiircd during consecutive 16 h with exposure times 
equal or twice the Pspin to smear out effects of the rota- 
tional pulsation (Table 1). The SWP and LWP exposures 
sample the orbital period. 

The spectra have been re-processed at VILSPA us- 
ing the lUE NEWSIPS pipeline used for the lUE Final 
Archive which applies the SWET extraction method as 
well as the latest flux calibrations and close-out camera 
sensitivity corrections (Garhart et al. 1997). Line-by-line 
images have been inspected for spurious features which 
have been identified and removed. 

2.3. The optical photometry 

Optical photometry was conducted in the period Octo- 
ber 18 23 1995 at the SAAO 0.75 m telescope and UCT 
Photometer employing a Hamamatsu R93402 GaAs pho- 
tomultiplicr. BVRI (Cousins) photometry was carried out 
on the first three nights performing symmetric modules 
with integration times of 30 s or 20 s respectively for the 
B and I, or V and R filters. The typical time resolution 
for the sequence of all four filters was ~12Gs, with inter- 
ruption every ^ lOmin for sky measurements. The times 
and durations of individual runs are reported in Table 1. 
The orbital period has not been fully sampled. 

Additional fast photometry in white light has been car- 
ried out on October 21 and 23, 1995 using the same pho- 
tometer, but employing a second channel photomultiplier 
for monitoring of a nearby comparison star. Continuous 
5 s integrations were obtained, with occasional interrup- 
tions (every 15-20 min) for sky measurements (lasting ~ 
20-30 s). During the two nights, the observations were car- 
ried out for 4.3 h and 2.5 h respectively. 

The photometric data have been reduced in a stan- 
dard manner with sky subtraction, extinction correction 
and transformation to the standard Cousins system using 
observations of E-region standards obtained on the same 
night. 



3. The UV spectrum of FO Aqr 

The grand average UV spectra of FO Aqr as observed with 
FOS and lUE are shown in Fig. 1 (upper panel), where the 
above flux difference is apparent. 

The UV hmiinosity in the 1150 -3200 A lUE range is 
6 X 10"^^ ergs s^^assuming a distance of 325 pc (Paper 1), a 
factor 1.3 larger than during previous UV observations in 
1990. The optical photometry also indicates a brightening 
of 0.17mag between the two epochs, indicating long term 
luminosity variations (see also sect. 8.5) 

The UV spectrum of FOAqr is typical of magnetic 
CVs (Chiappetti et al., 1989; de Martino 1995) with strong 
emissions of NV A1240, Si IV A1397, CIV A1550, HeH 
AA1640,2733 and Mgll A2800. The weaker Si IV with re- 
spect to NV emission, classes the IP nature (de Martino 

1995) . Weaker emissions from lower ionization states of 
different species such as CIII A1176, the blend of Si III 
A1298 multiplet, SiIIA1304 and geocoronal OIA1305, 
Silll A1895,SiH A1808, NIV A1718, NIII] A1747-1754, 
Aira A1855andAlII A1670, as weU as Hell A2307, possi- 
bly blended with C III A2297 , and Hell A2386, are identi- 
fied in the higher quality FOS spectrum. Also weak oxygen 
lines of OIV A1343 and O V A1371 are detected. Some of 
these lines are also observed in the HST/FOS spectra of 
AEAqr (Eracleous & Horne 1994), DQHer (Silber et al. 

1996) and PQGem (Stavroyiannopoulos et al. 1997). The 
presence of high ionization species together with extremely 
weak emissions (E.W. < 1 A ) of lower ionization species 
are characteristic of a higher ionization efficiency in IPs 
with respect to Polars (de Martino 1998). The line ratios 
NV/SilV and NV/Hell, when compared with photoion- 
ization models developed by Mauche et al. (1997) are close 
to the predicted values for an ionizing blackbody spectrum 
at 30 eV. 

In contrast to the lUE spectra, the FOS data allow 
us to finally detect the intrinsic Lya A1216 line. This ap- 
pears to be composed of a relatively deep (E.W.=5.4±0.1 
A) absorption and a weak emission (E.W.=1.4±0.lA ). 
The center wavelength of the absorption feature is how- 
ever red-shifted by ~4A with respect to rest wavelength 
and other emission line positions, while the weak emission 
is blue-shifted at 1206 A in the grand average spectrum. 
Discussion on the nature of this feature is left until sects. 5 
and 7, however the Lya absorption provides an upper limit 
to the hydrogen column density along the line of sight to 
FO Aqr. 

A pure damping Lorentzian profile (Bohlin 1975) con- 
volved with a 7A FWHM Gaussian has then been fit- 
ted to the Lya absorption line (Fig. 1, bottom panel). 
The resulting neutral hydrogen column density is iVu = 
(5.0 ± 1.5) X lO^^cm"^. The residual from the fit shows an 
emission line with maximum flux at ~ 1215 A , probably 
geocoronal or intrinsic, with an excess of flux in the blue 
wing possibly due to emission of Si III A 1206. 
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1500 2000 
Wavelength (A) 



2500 




Wavelength (A) 

Fig. 1. Upper panel: Grand average FOS (solid line) 
and lUE (dotted line) spectra of FOAqr in the range 
1150-2500A together with line identification. Bottom 
panel: Enlargement of the Lya feature in the average 
FOS spectrum. The solid line represents the observed 
spectrum, the dotted line a damped Lorentzian profile 
convolved with a 7A FWHM gaussian corresponding to 
A^jj = 5 X 10^°cm~^ and the dashed line represents the 
ratio between the two, enhancing the emission feature at 
1215A . 



The derived value for A^h is consistent with the to- 
tal interstellar column density in the direction of FO Aqr 
as derived from Dickey & Lockman (1990) and with the 
upper limit estimated from X-rays (Mukai et al. 1994). 
Assuming an average gas-to-dust ratio (Shull & van Steen- 
berg 1985), this upper Hmit corresponds to a reddening of 
Eb-v ~ 0.1. Although FOAqr was already known to be 
negligibly reddened from lUE observations (Chiappetti et 
al. 1989), an upper limit Eb-v = 0.013± 0.005 is derived 
from the absence of the 2200 A absorption in the FOS 
data. This indicates that, despite the coincidence, most of 
the neutral absorption is tmrclatcd to the interstellar dust 
and hence it is likely located within the binary system. 



4. Time series analysis 

The presence of periodicities in FO Aqr has been investi- 
gated in the FOS continua, emission lines and zero-order 
light as well as in the optical photometric data. 



4.1. HST UVdata 

Fluxes in five line-free continuum bands have been mea- 
sured in each FOS spectrum in the ranges A1265-1275, 
A1425-1450, A1675-1710, A2020-2100, A2410-2500. Line 
fluxes of Hell A1640, NV, Si IV and CIV and Lya have 
been computed adopting a method which uses for the 
continuum a power law distribution as found from a fit 
in the above continuum bands. Furthermore, since the 
low spectral resolution of FOS data prevents the study 
of UV line profiles, measures of the V/R ratios of emis- 
sion lines have been used to investigate possible mo- 
tions in the lines. These are defined as the ratios be- 
tween the integrated fluxes in the violet and red por- 
tions of the emission lines assuming as centroid wavelength 
that measured in the average profile. Such analysis is re- 
stricted to the strong emissions of Hell, CIV and NV 
lines whose FWZI are ±3000kms-\ ±4000kms-iand 
±3000 km s~^ respectively. 

In order to detect the active frequencies in the power 
spectrum and to optimize the S/N ratio, a Fourier analysis 
has been performed using the DFT algorithm of Deem- 
ing (1975) on the total UV continuum (sum of the five 
bands) and line (sum of NV, Si IV, CIV and Hell) fluxes 
and zero-order light. From Fig. 2, the dominance of the fl 
variability is apparent, being about twice the oj signal, as 
well as the presence of substantial power at the sideband 
and orbital harmonic frequencies. To distinguish real sig- 
nals from artifacts due to the sampling of the HST orbit, 
a least-square technique was applied to the each data set 
which fits simultaneously multiple sinusoids at fixed fre- 
quencies. A synthetic light curve with the same temporal 
sampling of the data was created and subtracted (resid- 
uals). The continuum and zero-order light reveal, besides 
the uj and frequencies, also the a; — 2 f7, lu — CI and fl + uj 
sidebands, whereas in the emission lines the 2 Q and co — Q 
frequencies are detected. In Fig. 2 the amplitude spectra 
relative to the multiple sinusoids and the residuals arc also 
shown for comparison. It should be noted that peaks at 
~ 54day~^ and ~ 84day~^ are identified as sideband fe- 
qtiencies of the spin and HST orbital frequencies. These 
are removed by the method as shown by the residuals. 
Then a five (four) frequency composite sinusoidal func- 
tion for each spectral band (emission line) has been used 
and the derived amplitudes, reported in Table 2, have been 
compared with the average power in the DFTs of the resid- 
uals (cr) in the range of frequencies of interest (i.e. v ^ 
1.4 mHz), (column 8). While the signals at w — 2f2 (con- 
tinuum) and at 20 (emission lines) on average fulfil a 4 
a criterium, the other sidebands are between between 2.1 
and 3.5 a. 

A further check has been performed using the CLEAN 
algorithm (Roberts et al. 1987) which removes the win- 
dowing effects of the HST orbit. The CLEANED power 
spectra, adopting a gain of 0.1 and 500 iterations, indeed 
reveal the presence of the iv — 20, and iv — fl sidebands 
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Fig. 2. From top to bottom: DFTs of the UV contimmm, zero-order light, emission hne fluxes and V/R ratios of Hell, 
C II and N V (solid line) . DFTs of the sinusoidal funetions (dotted line) together with those of the residuals (dashed 
line) are also shown. The frequencies used in the multiple sinusoidal fits are marked with vertical solid lines. 



in the continuum and the 211 in the emission lines. The 
lack of significant power at the other frequencies is consis- 
tent with the previous analysis. Hence, these weakly active 
frequencies will be considered with some caution in this 
analysis. 

A strong colour effect in the UV continuum is detected 
with amplitudes decreasing at longer wavelengths. Differ- 
ent from other lines is the behaviour in the Lya feature, 
whose absorption component is modulated at the fl fre- 
quency, while a variability at the spin is at a 2a level. No 
significant variations are detected in the equivalent width 
of the absorption as well as in the emission component. 

In contrast to the flux behaviour, the V/R ratios are 
variable only at the u) frequency (Fig. 2, bottom panels). 



Noteworthy is the marginal spin variability in the CIV 
line. The amplitudes of the spin modulation, obtained 
from the least-square fits, are reported in the last column 
of Table 2. 



4-2. Optical data 

The analysis of BVRI data acquired during three nights 
has been performed following the same procedure adopted 
for the HST data. Contrary to previous observations, the 
orbital variability also dominates in the optical being ~ 1.5 
times the spin modulation. The presence of sideband mod- 
ulations is more uncertain because of the lower quality of 
the data. Nevertheless, to ensure uniformity between UV 
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Table 2. Amplitudes of the modulations for the UV continuum, optical and UV emission lines derived from a multi- 
component sinusoidal fit. 



Band (A) 


An 


Am 


Au; 


A(^_2n 


Aoj-n 


Aoj+n 


Ares A^{yiR) 


1265 — 1275 


1.730(52) 




0.893(51) 


0.913(72) 


0.391(55) 


0.693(72) 


0.181 


1425 - 1450 


1.784(42) 




0.867(43) 


0.694(59) 


0.453(44) 


0.441(61) 


0.166 


1675 1710 


1.779(37) 




0.738(37) 


0.685(51) 


0.344(37) 


0.420(52) 


0.157 


2020-2100 


1.416(29) 




0.656(29) 


0.553(40) 


0.303(30) 


0.334(41) 


0.125 


2410-2500 


1.130(25) 




0.627(25) 


0.463(34) 


0.236(25) 


0.257(34) 


0.108 


2900-2985 


0.942(45) 














Zero-Order 


0.578(15) 




0.417(14) 


0.323(21) 


0.120(15) 


0.178(21) 


0.062 


B 


0.354(25) 




0.259(38) 


0.025(33) 


0.046(35) 


0.047(31) 


0.045 


V 


0.182(14) 




0.129(22) 


0.050(19) 


0.006(21) 


0.023(18) 


0.026 


Rc 


0.134(13) 




0.091(20) 


0.026(16) 


0.020(21) 


0.037(16) 


0.021 


Ic 


0.109(10) 




0.054(16) 


0.022(14) 


0.031(17) 


0.015(13) 


0.016 



NVflux 14.37(1.98) 6.75(1.97) 12.75(1.37) 12.25(1.37) 3.32 0.256(41) 

Si IV flux 15.40(1.54) 13.17(1.53) 7.60(1.07) 7.77(1.07) 2.75 

CIV flux 77.10(2.87) 36.74(2.86) 35.71(1.98) 22.22(1.99) 16.99 0.057(08) 

Hell flux 31.81(1.59) 14.68(1.60) 22.57(1.10) 13.40(1.11) 2.75 0.173(10) 

Lya abs. 23.35(1.08) 6.83(1.07) 3.41 



Notes: (1) Continuum flux amplitudes are in units of 10~^* ergs cm~^s~^ A ^while line fluxes are in units of lO"^'' ergs cm~^s 

Errors in parentheses arc referred to the last significant digits. 

(2) Ar-es is the average amplitude of the DFT of the residuals, calculated for frequencies ^1.4 mHz. 

(3) For the AA 2900-2985 band only the orbital amplitude is reported as derived from lUE data. 



and optical results, a least-square fit to the data has been 
applied using the same five frequency sinusoidal function, 
i.e. n, w, w — 2 n, w — f2, w + f2. The resulting amplitudes, 
reported in Table 2, when compared with the noise in the 
residuals (column 8), can be considered as upper limits. 

As far as the fast photometry is concerned, the low 
quality of the data only allows the detection of the spin 
and orbital variabilities. In particular, the latter is de- 
tected on the first night with a pronounced dip which is 
not consistent with the refined orbital ephemeris based on 
orbital minima recently given by Patterson et al. (1998), 
which defines the inferior conjunction of the secondary 
star. Therefore these data will not be used for a multi- 
wavelength analysis of the pulsations. 

/^From this analysis new times of maxima for the or- 
bital and rotational modulations are derived for the UV 
continuum and optical light: 

iIJ£>X = 2 449 971.2251 ± 0.0006 in the UV; 

ilJD^t? = 2 449 971.11181 ± 0.00010 in the UV; 



J£>™|^ = 2 449 971.2235 ± 0.0034 in the optical; 
if JL>^t^ = 2 449 971.11021 ± 0.00035 in the optical; 

Both UV and optical orbital maxima lead by A$orb 
= 0.145 those predicted by Patterson et al.'s ephemeris. 
Such phase difference, discussed in sect. 7, is consistent 
with the previous UV results (Paper 1). 

On the other hand, the optical rotational maximum 
agrees within 8 per cent with that predicted by the new 
revised cubic ephemeris given by Patterson (1998, private 
communication) : 

HJD 2 444 782.9168(2) + 0.014519035(2) 

7.002(7) 10-^3 E'^ - 1.556(2) lO'^^ E'^ (1) 

The UV rotational pulses lag the optical by A$rot = 
0.186. Such colour effect will be discussed in more detail 
in sect. 5. 

Furthermore, the time of coherence between spin and 
beat modulations is found to be HJDr.^ = 2449971.2335 
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± 0.0070 in both UV and optical. Throughout this pa- 
per the Patterson ct al.'s orbital ephemeris will be used 
but $orb = 0.0 will refer to the orbital maximum. Hence, 
phase coherence occurs at $orb = 0.90±0.03. i.e. close to 
the orbital maximum, whilst in 1988 and 1990 it was found 
close to the orbital minimum (Osborne & Mukai 1989; Pa- 
per 1). Such phase changes are not uncommon for FO Aqr 
(Semeniuk & Kaluzny 1988; Hellier et al. 1990). 

5. The rotational modulation 

UV continuum, zero-order light, and UV emission line 
fluxes as well as their V/R ratios have been folded in 56 
phase bins along Pspin- The spin pulses, pre-whitened from 
all other frequency variations, are shown in Fig. 3 together 
with the optical B band pulses folded in 28 phase bins. A 
strong colour effect is observed in both amplitudes and 
phasing. Fractional amplitudes (amplitudes of the sinu- 
soid Asin(a;t -|- </>) divided by the average value) decrease 
from 26% in the far-UV to 16% in the near-UV and to 
10% in the optical. The far-UV pulse maximum is broader 
and lags by A(/)spi„ = 0.091 ± 0.010 and 0.186±0.008 the 
near-UV and optical maxima respectively. While the UV 
line fluxes follow the near-UV modulation, their V/R ra- 
tios are in phase with the far-UV, with a maximum blue 
shift when the far-UV pulse is at maximum. The pulsa- 
tion in the line fluxes are 18%, 12%, and 10% in He II, 
CIV and Si IV. The V/R ratios are generally <1 indi- 
cating the presence of a dominant blue component in the 
lines, which is also visible from the extended blue wings in 
the average spectrum. The velocity displacements indicate 
the presence of a spin S-wave in the line profiles similar 
to the optical (Hellier et al. 1990). Both continuum and 
emission lines therefore strongly indicate the presence of 
two components which affect the rotational modulation in 
FO Aqr. 

The 797 FOS spectra have been spin-folded into 20 phase 
bins. A total of 780 light curves, each sampling a wave- 
length bin of 1.8 A , were then produced and fitted 
with a sinusoid. The resulting amplitudes define the ro- 
tational pulsed spectrum as Fa = 2Aa. This spectrum 
(shown in the enlargement of the lower panel of Fig. 3), 
has to be regarded as an upper limit to the modulated 
flux since no pre whitening could be performed given the 
low S/N of each 780 light curves. This spectrum gives 
evidence of modulation not only in the main emission 
lines and Lya absorption but also in the weaker emis- 
sion features identified in sect. 3. Broad band UV con- 
tinuum and optical photometric spin pulsed fiuxes, ob- 
tained from the multi-frequency fit and reported in Fig. 3, 
provide a correct description of the rotational pulsed en- 
ergy distribution. A spectral fit to the broad band UV 
and optical spectrum, using a composite spectral func- 
tion, consisting of two blackbodies, gives 37500±500K 
and 12 000±400K (x^d = 0.91). The projected fractional 
area of the hot component is ~0.11 A^j, while the cool one 



covers ^ 10.7 A^d, for R^d = 8 x 10® cm, and d=325pc 
(Paper 1). 

The FOS rotational pulsed spectrum shows the pre- 
sence of a Lya absorption feature which gives a hydrogen 
column density of 8 ± 2 x 10^" cm^^. On the other hand, 
assuming A^h = 5 x 10^" cm^^, as derived from the grand 
average spectrum, a composite function consisting of a 
white dwarf model atmosphere with at 36 000 K and of 
the same 12 000 K blackbody gives an equally satisfactory 
fit to the whole FOS spectrum. For a distance of 325 pc 
the radius of the white dwarf is 4.9x10^ cm, in agreement 
with that of the hot blackbody component. 
While the detection of the hot component is new, a com- 
parison with previous spectral analysis of the optical and 
IR spin pulses observed in 1990 shows that the tempera- 
ture of the cool component has not changed with time but 
instead it suffered a decrease in area by a factor of ~ 1.5. 

6. The sideband modulations 

The beat lu — Q modulation, although weak among the 
detected sidebands, shows an anti- phased behaviour be- 
tween line fluxes and continuum (Fig. 4a). The UV line 
fluxes, with average fractional amplitudes of ^ 11%, show 
a minimum when the UV continuum is maximum. Their 
maximum shows a dip-like feature centered on the UV 
continumm minimum. Colour effects are also encountered, 
the far-UV pulses being stronger (^ 11%) than the near- 
UV (^ 6%) ones and lagging by ^ 0.2 in phase the near- 
UV. The trend of a decrease in the amplitudes at increas- 
ing wavelengths is confirmed in the optical where an upper 
limit of ~ 2% can be set to the fractional amplitudes. 
A similar colour effect is also observed in the lu + Cl side- 
band pulsation with fractional amplitudes ranging from 
~ 20% in the far-UV to ^ 7% in the near-UV, displaying 
a broadening of the maximum towards longer wavelengths 
(Fig. 4b). 

The strong uj — 2Q pulsation only shows a wavelength de- 
pendence in the fractional amplitudes, being in the far-UV 
^ 26% and decreasing to ^ 12% in the near-UV (Fig. 4c) 
and to ^ 2 — 3% (upper limit) in the optical. 
The modulated spectra at the three frequencies have been 
derived following the same procedure as for the rotational 
pulsed spectrum. The anti-phased behaviour of the emis- 
sion lines is seen in the beat pulsed spectrum, where these 
are seen as absorption features, except for residuals in the 
CIV and Hell fines (Fig4d, bottom panel). The modula- 
tion spectra at lu — 2fl and lu + frequencies show weak 
emissions at Si IV, CIV and Hell indicating a marginal 
variability at these frequencies. The UV continuum energy 
distributions of these variabilities are best represented by 
power laws Fa oc A^" with spectral index a ~ 0.8 — 1.7, 
rather than blackbodies (20000 - 25 000 K), possibly sug- 
gesting that more than one component is acting. Given 
the low level of confidence of the sideband variabilities it 
is not possible to derive further information. 
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Fig. 3. Upper left: Spin pulses in the far-UV, near-UV, zero order and B band (upper left panel). Upper right: 
CIV and Hell flux and V/R rotational curves. Fluxes are fractional as described in the text, the average value has 
been subtracted. Bottom panel: Rotational broad band UV and optical energy distribution together with the best fit 
composite function: a hot (37 500 K) (dotted line) and a cool (12 000 K) (long-dashed line) blackbody function (sum: 
solid line). A composite function (solid line) consisting of a 36 000 K white dwarf model spectrum and of the same 
12 000 K blackbody absorbed by iVn = 5 x 10^°cm~^ is shown in the insert figure together with the FOS rotational 
pulsed spectrum described in the text. 



7. The orbital vEtriability 

The UV and optical orbital modulations have been in- 
vestigated folding the FOS continuum broad band and 
emission line fluxes in 28 orbital phase bins. The light 
curves have been prewhitened by the other active frequen- 
cies using the results of the multi-frequency fits. For the 
lUE data, continuum broad band and emission line flux 



measures have been performed on each SWP and LWP 
spectrum. Three broad bands have been selected in each 
spectral range, five of them coinciding with the FOS se- 
lected bands and a sixth one in the range A 2900-2985. 
The contribution of the spin pulsation, as derived from 
the multi-frequency fit has been removed. The best fit 
blackbody spin pulsed spectrum has been used to allow 
prewhitening in the range A2900-2985. The FOS and lUE 
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Fig. 4. a The beat w — ^l modulation in the far-UV, near-UV continua and line fluxes of C IV and He II. b The a; + f2 
and c the ui — 2Sl light curves in three UV bands and zero order light. Fluxes are fractional, d From bottom to top, 
the pulsation spectra at the beat w — ^i, w — 2 and u) + Cl frequencies (solid line) together with broad band UV fluxes 
derived from the multi- frequency fits (filled circles). 



broad band continuum fluxes in the far-UV, mid-UV and 
near-UV as well as the zero order and B band light curves 
are reported in the left panel of Fig. 5, while the emis- 
sion line fluxes of Si IV, C IV and He II arc shown in the 
right panel. The orbital gaps due to the HST sampling are 
apparent. 



A strong colour dependence is encountered in the mod- 
ulation amplitudes as well as the phasing. Fractional am- 
plitudes range from 40% in the far-UV to 28% in the near- 
UV (lUE band) and 15% in the optical. The modiflation 
amplitudes then have increased by a factor ~ 2.5 in the UV 
and ~ 1.5 in the optical with respect to 1990. However, the 
phasing of UV maximum and minimum has not changed 
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Fig. 5. Upper left: Orbital modulation in the far-UV, mid-UV and near-UV, zero order and B light (upper left panel). 
Upper right: Si IV, CIV and Hell flux light curves. Fluxes are fractional as described in the text and the average 
value has been subtracted. The lUE measures are reported together with their phase coverage. Bottom panel: The 
orbital broad band UV and optical modulated fluxes represented with the best fit (solid line) composite function as 
described in the text. The two blackbodies are represented with dotted lines. A hot (21500 K) blackbody function 
describes the UV FOS spectrum (shown in the inserted figure.) 



with time, occurring at $orb = 0.86 and $orb = 0.34, re- 
spectively. The UV modulation is more sinusoidal, whilst 
the optical light curve is more structured with a double 
humped maximum between $orb = 0.75 and $orb = 0.0. 
A comparison with the optical behaviour in 1990 indi- 
cates an absence of a broad maximum centered at 'I'orb = 
0.0 and a less defined minimum. The current observations 
are inadequate to resolve the orbital dip due to a graz- 



ing eclipse of the accretion disc in either UV and optical 
ranges. 

The orbital modulation in the UV emission line fluxes is 
strong with fractional amplitudes of 40% in NV, 20% in 
Si IV, 29% in CIV and 23% in Hell and almost in phase 
with the UV continuum. 

The spectrum of the UV orbital variability derived 
with the same procedure as described before is shown in 
the enlargement of Fig. 5 (bottom panel). ^From the in- 
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serted figure the Lya absorption feature is apparent and 
is consistent witlr the neutral hydrogen column density of 
8 ± 2 X 10^° cni^^ inferred from the spin pulsed spectrum. 
Hence, while this absorption in the orbital pulsation spec- 
trum is clearly circumstellar, the same nature in the spin 
pulsed spectrum cannot be excluded. 

The UV FOS spectrum requires a hot component with 
a blackbody temperature of 21 500±500K (Fig. 5 bottom 
panel enlargement). The composite UV and optical broad 
band energy distrubution confirms the previous results on 
the presence of two components, a hot at 19 500 ± 500 K 
and a cool one at 5 700 ± 200 K (x^^d = 1-1) (Fig- 5, bot- 
tom panel). With current UV observations, it is now pos- 
sible to constrain the temperature of the hot emitting re- 
gion. The temperature of the cool component is in agree- 
ment, within errors, with that inferred in Paper 1. A sub- 
stantial increase by a factor of ~ 2.6 in the area of the hot 
region is found when compared to the 1990 epoch, which 
is 12 Awd- The emitting area of the cool component is 
instead similar to that previously derived (Paper 1). 

8. Discussion 

The HST/FOS and lUE spectroscopy has revealed new 
insights in the UV variability of FO Aqr. 

8.1. The periodic variations 

The UV continuum and emission line fluxes are found to 
be strongly variable at the orbital period. This periodic- 
ity also dominates the optical range where FO Aqr was 
previously found to be spin dominated. The time series 
analysis indicates the presence of other periodicities, the 
negative lu — 2fl sideband being much stronger than the 
beat w — f7 in both UV and optical ranges. The presence 
of sidebands with different amplitudes at different epochs 
is rather common in FO Aqr, although the beat is usually 
the strongest (Patterson & Steiner 1983; Warner 1986; Se- 
meniuk & Kaluzny 1988, Chiappetti et al. 1989; Paper 1, 
Marsh & Duck 1996). In particular the intermittent oc- 
currence of a stronger pulsation at w — 2 f2 frequency was 
already noticed (Warner 1986; Patterson et al. 1998). The 
strong negative sideband u — 2Q, cannot be produced by 
an amplitude modulation at 2il frequency of the rotational 
pulses, since the positive sideband co + 2^1 should have 
been present. Also, an orbital variability of the amplitude 
of the w — n modulation cannot be responsible alone since 
it is too weak. Hence the u} — 2fl pulsation should be dom- 
inated by the effects of an unmodulated illumination from 
the white dwarf, which naturally gives rise to the orbital 
variability (Warner 1986). 

The occurrence of coherence between spin and beat pulsa- 
tions appears to be different from epoch to epoch (Seme- 
niuk & Kaluzny 1988; Osborne & Mukai 1989; Paper 1). 
It was proposed that phase coherence close to the opti- 
cal orbital minimum could be possible if the reprocessing 



site(s) are viewing the lower accreting pole (Paper 1). The 
observed shift of half an orbital cycle would then imply 
that the reprocessing region(s), are now viewing the main 
accreting pole, as predicted by the standard reprocessing 
scenario (Warner 1986). 

The behaviour of UV emission lines is different between 
fluxes and V/R ratios. The line fluxes are strongly vari- 
able at the orbital period, the spin variability being 1.6 
times lower. On the other hand, their V/R ratios only 
show a rotational S-wave, but that of C IV line is surpris- 
ingly weak. The lack of detection of an orbital S-wave, 
can be ascribed to the low amplitude (^ 300 400 kms^^) 
velocity displacements known from optical data (Hellier 
et al. 1989; Marsh & Duck 1996), which are not detected 
because of the low spectral resolution of the FOS data. 

8.2. The rotational pulses 

Both shapes and amplitudes of UV and optical contin- 
num spin pulses indicate the presence of two components, 
one dominating the near-UV and optical ranges, already 
identified in Paper 1 and a new contribution dominating 
the far-UV pulses which lags by ^ 0.2 in phase the first 
one. Furthermore a different behaviour between emission 
line fluxes and V/R ratios is observed. While the latter 
show a spin S-wave in phase with the far-UV continuum, 
the line fluxes follow the near-UV and optical pulsations. 
The maximum blue-shift found at rotational maximum 
of the far-UV pulses indicates that the bulk of velocity 
motions in the emission lines maps the innermost regions 
of the accretion curtain. The outer curtain regions are 
then responsible for X-ray illumination effects seen in the 
line fluxes and near-UV and optical continua. A direct 
comparison with previous X-ray observations reported by 
Beardmorc et al. (1998) is not possible since, adopting 
their linear spin ephemeris, the UV and optical maxima 
lag by A(E>spin=0.4 their predicted optical maximum. How- 
ever, the X-ray pulse maxima observed by Beardmore et 
al. (1998), typically lag by A$spi„ =0.2 their optical phase 
zero (see their Fig. 3), consistently with the lag observed 
between the far-UV and optical pulses. Hence this differ- 
ence is an indication that the far sides of the accretion 
curtain come into view earlier than the innermost regions. 

The spectrum of the pulsation reveals regions at ~ 
37 000 K covering a relatively large area, ^ 0.1 A„d, with 
respect to typical X-ray fractional areas / < lO^"* (Rosen 
1992). Such hot components have been also observed in 
the IPs PQ Gem (Stavroyiannopoulos et al. 1997) and 
EXHya (de Martino 1998). The presence of the Lya ab- 
sorption feature in this spectrum can be partially due to 
the photospheric absorption of the heated white dwarf 
with similar temperature and fractional area as a black- 
body representation. On the other hand, both orbital and 
rotational modulated spectra give similar values of A^h if 
this absorption is of circumstellar nature. Hence, while 
only in AE Aqr the UV pulses are clearly associated with 
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the heated white dwarf (Eracleous & Home 1994, 1996), 
those in FO Aqr can be associated with either the inner- 
most regions of the accretion curtain onto the white dwarf 
or its heated polar regions. 

The second component identified as a cool 12 000 K re- 
gion covers ^ llAwd, a factor ~ 1.5 lower than previously 
found in Paper 1. Thus the decrease in the optical ampli- 
tudes does not involve substantial changes in temperatures 
but in the size of the accretion curtain. Such lowc;r temper- 
atures characterizing the near-UV and optical pulses are 
also recognized in other IPs (de Martino et al. 1995; Welsh 
& Martell 1996; Stavroyiannopoulos et al. 1997; de Mar- 
tino 1998). Although a two component pulsed emission 
might be a crude representation, it is clear that temper- 
ature gradients are present within the accretion curtain 
extending up to ~ 6Rwd- 

The bolometric flux involved in the spin mod- 
ulation due to both components amounts to 
6.3xl0~^^ ergs cm~ ^ s~ ""^ . Although no contemporary 
X-ray observations are available this accounts for ~ 26% 
the total accretion luminosity as derived from ASCA 
1993 observations (Mukai et al. 1994). 

8.3. The sidebands variability 

The UV continuum pulsations observed at the sideband 
frequencies, uj — 217, u) — VI and w -I- fi, indicate the pres- 
ence of a relatively hot component ~ 20 000-25 000 K. The 
lack of adequate data in the optical range does not allow 
one to confirm the cool 7000K), and hence possible 
second component, in the beat pulsed energy distribution 
as found in Paper 1. The phase lags of the far-UV max- 
imum with respect to that in the near-UV in the iv — Q. 
and (jj -\-^ pulsations are similar to that observed in the 
rotational pulses. This is consistent with the pulsations 
being produced by amplitude variations of the spin pulses 
at the orbital period. In contrast, the prominent negative 
sideband u) — 20. variability is not affected by phase shift 
effects, indicating that indeed such variability is mainly 
due to an aspect dependence of the reprocessing site at 
the orbital period. 

No strong pulsation in the UV emission lines is observed 
at these frequencies except for the interesting anti-phased 
behaviour of these lines at the beat period. These are ob- 
served as weak absorption features in the modulated spec- 
trum. Such behaviour, although much more prominent, is 
also observed in the IP PQGem (Stavroyiannopoulos et 
al. 1997). Though this is not easy to understand, a pos- 
sibility could be that the reprocessing site, producing the 
LO — Vl component in the emission lines, is viewing the lower 
pole instead of the main X-ray illuminating pole as also 
suggested by Stavroyiannopoulos et al. (1997). 



8-4-. The orbital variability 

The present study confirms previous results where the or- 
bital modulation is composed by two contributions, iden- 
tified as the illuminated bulge and the heated face of the 
secondary star, the former being at superior conjunction 
at (E>orb = 0.86, while the latter is at superior conjunc- 
tion at 3>orb = 0.0. The double-humped maximum in the 
optical light curve can be understood in terms of relative 
proportion of a strong bulge contribution with respect to 
that of the secondary star. Indeed no changes in the tem- 
peratures are found (the hot one is better constrained with 
the present data), but a substantial change by a factor of 
~ 2.6 since 1990 is found in the emitting area of the bulge 
itself. All this indicates that illumination effects are ba- 
sically unchanged whilst the inflated part of the disc has 
increased. 

The total bolometric flux involved in the orbital 
variability amounts to 2.7xl0~^''ergs cm~^s~^which is 
a factor ~ 4 larger than that of the rotational pul- 
sation. Neglecting the sideband contributions at first 
approximation, the total modulated flux amounts to 
3.3x10"^° ergs cm~^s~^and corresponds to a reprocessed 
luminosity of 4.2 x 10'^^ ergs s~^which is approximately the 
order of magnitude of the accretion luminosity derived 
from X-rays (Mukai et al. 1994). Then assuming the bal- 
ance of the energy budgets of the reprocessed and primary 
X-ray radiations, and hence of the accretion luminosity, an 
estimate of the accretion rate of M = 6.7 x 10~^° Mq yr~^ 
is derived. 

8.5. The long term variability 

FO Aqr has displayed a change in its power spectrum at 
optical and UV wavelengths on a time scale of five years. It 
was also brighter by 0.3 mag and 0.2 mag in the two ranges 
with respect to 1990. This difference is accounted for by 
the orbital modulated flux. The study of the spectra of the 
periodic variabilities has shown that a shrinking of the ac- 
cretion curtain by a factor of 1.5 has occurred while the 
inflated part of the disc has increasc^d in area by a factor of 

2.6. Such changes indicate variations in the accretion pa- 
rameters. Worth noticing is that the unmodulated UV and 
optical continuTim component has not changed with time, 
indicating that the steady emission from the accretion disc 
(Paper 1) has not been affected. These results are in agree- 
ment with the long term trend of the X-ray power spectra 
(Beardmore et al. 1998) which showed that FOAqr was 
dominated by the spin pulsation in 1990, while in 1988 
and in 1993 prominent orbital and sideband variabilities 
were present. Changes from a predominant disc- fed ac- 
cretion to a disc-overflow (or stream-fed) accretion have 
been the natural explanation for such changes in the X-ray 
power spectra. As proposed in Paper 1, the bulge provides 
a source for a variable mass transfer onto the white dwarf. 



de Martino et al.: Time Resolved HST and lUE 



13 



and an increase in its dimensions accounts for a predomi- 
nant disc-overflow towards the white dwarf. 

Beardmore et al. (1998) suggested that changes in the 
accretion mode could be triggered by variations in the 
mass accretion rate and the analysis presented here is in- 
deed in favour of this hypothesis. An estimate of changes 
in the accretion rate producing a shrinking of the accre- 
tion curtain can be inferred by the relation betwen mag- 
netospheric radius, accretion rate and magnetic moment 
(Norton & Watson 1989): 

Tmas = </-2.7 X IQiV^/'Mre'/'Mj]/' 

(2) 

where < 1 is a dimensionless factor accounting for 
a departure from spherical symmetry, ix^^ is the magnetic 
moment in units of 10^^ Gcm^, Mie is the mass accretion 
rate in units of 10"'^^ gs"-'^ and M^d is the white dwarf 
mass in units of M© . Hence assuming that the accretion 
curtain reaches the magnetospheric boundary, a reduction 
of the linear extension by a factor of ^ 1.2 implies that 
the accretion rate has increased by a factor ~ 2 in five 
years. 

The long term spin period variations observed from 1981 
to 1997 in FO Aqr (Patterson et al. 1998), which changed 
from a spin-down to a recent spin-up since 1992, were pro- 
posed to be due to variations around the equilibrium pe- 
riod produced by long term variations in M. The increase 
in brightness level and the results found in the present 
analysis are strongly in favour of this interpretation. 

9. Conclusions 

The present HST/FOS and lUE spectroscopy has allowed 
us for the first time to infer the characteristics of multiple 
periodicities in the UV emission of FO Aqr. Also, coordi- 
nated optical photometry has provided further constraints 
summarized as follows: 

1) The rotational pulsations in the UV and optical 
are consistent with the current accretion curtain sce- 
nario. They reveal the presence of strong temperature 
gradients within the curtain moving from the footprints 
onto the white dwarf surface to several white dwarf 
radii. A spin S-wave is observed in the UV emission lines 
which maps the innermost regions of the accretion curtain. 

2) Orbital sideband variabilities indicate that repro- 
cessing is occurring at fixed regions within the binary 
frame. Different behaviour is observed in the emission 
lines and continuum beat pulses indicating a complex 
reprocessing scenario. 

3) The orbital UV variability is confirmed to arise 
from the illuminated side of the inflated disc (bulge) while 
the optical modulation is produced by heating effects at 



the secondary star. 

^) A change in the relative proportion of rotational 
and orbital modulation amplitudes is found on a timescale 
of five years. These are interpreted as a reduction in the 
dimensions of the accretion curtain accompained by an 
increase in the bulge size. 

These observations indicate a long term change in the 
accretion mode where FOAqr has switched from a disc- 
fed to a disk-overflow state triggered by changes in the 
mass accretion rate. 
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